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Abstract 

This  report  develops  a  variational  technique  for  the  analysis  of 
the  vibration  characteristics  of  an  open  cylindrical  cantilevered  shell. 
The  technique  is  developed  by  modifying  Reissner's  principle,  which  nor¬ 
mally  applies  to  static  problems,  through  the  use  of  Hamilton's  princi¬ 
ple  so  that  it  applies  to  dynamic  problems.  The  variational  te-  .ique 
is  first  derived  in  general  for  an  elastic  system,  and  then  spet  j-cally 
tailored  to  an  open  cylindrical  cantilevered  shell.  The  techni 
implemented  by  first  finding  a  general  solution  which  satisfies  the 
equations  of  motion  for  a  cylindrical  shell.  A  method  is  then  formula¬ 
ted  to  use  this  general  solution  to  construct  a  set  of  trial  solution 
functions.  With  the  variational  method,  the  coefficients  to  this  trial 
solution  function  are  then  calculated  so  that  the  function  not  only 
satisfies  the  equations  of  motion,  but  also  the  boundary  conditions 
around  the  four  edges  of  the  shell.  A  computer  method  was  developed  to 
perform  the  necessary  calculations  to  implement  the  variational  pro¬ 
cedure,  but  preliminary  results  have  shown  that  numerical  problems  must 
be  eliminated  before  accurate  results  can  be  expected. 

Experimental  data  for  an  open  cylindrical  cantilevered  shell  was 
also  collected  on  a  modal  analyzer.  The  results  are  presented  and 
discussed. 


ix 


I.  Introduction 


Background 

The  vibration  characteristics  of  a  cantilevered  open  shallow  cylin¬ 
drical  shell  of  rectangular  planform  are  quite  different  than  those  of 
a  similarly  dimensioned  flat  cantilevered  plate.  The  determination  of 
these  differences  is  of  significant  technical  importance  in  areas  such 
as  in  the  analysis  of  turbomachinery  blades  where  accurate  natural  fre¬ 
quencies  and  mode  shapes  are  desirable.  Unfortunately,  tabulated  vibra¬ 
tion  data  is  not  as  readily  available  for  cantilevered  shells  as  it  is 
for  plates  and  beams.  In  fact,  in  Vibration  of  Shells  by  Leissa  (Ref  1) , 
one  of  the  most  comprehensive  collections  of  data  in  this  area,  there 
is  no  numerical  data  given  for  a  cantilevered  she^l  as  described  above. 
This  lack  of  data  is  a  result  of  the  equations  of  motion  becoming  more 
complex  in  going  from  a  plate  to  a  shell.  Thi3  increase  in  complexity 
allows  exact  solutions  to  but  a  few  special  cases  of  boundary  conditions . 
For  open  cylindrical  shells  some  of  these  special  cases  Include  those 
shells  which  have  shear  diaphragms  supporting  at  least  two  opposite 
edges.  Basically,  the  shear  diaphragm  is  the  analog  to  "simply  suppor¬ 
ted"  in  linear  beam  and  plate  theory.  A  shear  diaphragm  can  resist  all 
but  the  translation  normal  to  the  edge  and  the  bending  moment  along  the 
edge.  Therefore,  the  shear  diaphragm  does  not  represent  the  necessary 
boundary  conditions  needed  to  model  a  cantilevered  shell.  As  a  result, 
even  though  general  solutions  to  the  equations  of  motion  do  exist,  most 
practical  boundary  conditions  require  the  use  of  approximate  techniques 
such  as  the  Galerkin  or  the  Rayleigh-Ritz  methods.  Unfortunately,  these 
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methods  do  require  trial  functions  which  must  approximately  satisfy  at 
least  the  geometric  boundary  conditions.  This  requirement,  once  again 
restricts  the  class  of  solvable  problems  to  certain  types  of  boundary 
conditions.  To  overcome  this  restriction,  numerical  methods  such  as 
finite  elements  and  finite  differences  are  commonly  used.  In  an  effort 
to  develop  a  continuum  approach,  this  study  will  Investigate  a  varia¬ 
tional  method  derived  from  extending  Reissner's  principle  via  Hamilton's 
principle.  The  method  generates  a  solution  satisfying  the  boundary  con¬ 
ditions  by  a  superposition  of  a  set  of  functions,  each  of  which  satis¬ 
fies  the  governing  equations  of  motion,  but  not  necessarily  the  boundary 
conditions.  This  method  has  been  successfully  used  to  determine  the 
natural  frequencies  of  a  partially  clamped  circular  plate  (Ref  2) . 


Objective  and  Approach 

The  objective  of  this  study  is  to  determine  if  it  is  possible  and 
practical  to  apply  a  variational  procedure  to  calculate  the  free  vibra¬ 
tion  characteristics  of  shells  with  mixed  boundary  conditions.  In  par¬ 
ticular,  the  procedure  will  be  used  on  a  cantilevered  open  cylindrical 
shell . 

The  method  will  be  developed  using  thin  elastic  shell  theory 
specialized  to  the  case  of  cylindrical  shells,  and  then  applied  to  an 
open  cantilevered  cylindrical  shell.  A  procedure  will  then  be  devised 
to  calculate  the  coefficients  necessary  to  superimpose  trial  functions 
to  construct  an  approximate  solution.  A  computer  program  will  be  deve¬ 
loped  that  implements  the  variational  procedure.  Experimental  data  col¬ 
lected  on  a  modal  analyzer  will  also  be  presented  and  discussed. 


II.  Extension  of  Reissner's  Principle 


Reissner's  Principle 

Reissner's  principle  belongs  to  a  class  of  variational  techniques 
known  as  stationary  principles.  Unlike  the  principles  of  minimum  total 
potential  energy  and  minimum  total  complimentary  energy,  which  ensure 
the  extremum  found  is  a  minimum,  Reissner's  principle  ensures  only  that 
a  stationary  value  is  found.  The  exact  nature  of  the  extremum  is  not 
known.  Reissner's  principle  is  a  generalization  of  the  principle  of 
stationary  potential  energy  through  the  use  of  Lagrange  multipliers. 

The  principle  of  stationary  potential  energy  is  derived  from  the 
principle  of  virtual  work  by  assuming  that  the  applied  loads  are  con¬ 
servative  and  do  not  vary  in  magnitude  or  direction  during  virtual  dis¬ 
placements  (Ref  8:67).  This  principle  can  be  stated  as 


61-0 


(1) 


where 


Fiui 


Ujds 


(2) 


and 

WCu^)  -  Strain  energy  density  of  system  in  terms  of  dis¬ 
placements  and  their  derivatives 
u^  -  Displacements 

-  Body  forces  per  unit  volume 
T  *  *  Prescribed  surface  tractions 


V  ■  Total  volume  of  system 

*  Portion  of  boundary  where  tractions  are 
prescribed 

Reissner  generalized  this  functional  to  allow  simultaneous  artlbrary 
variations  in  displacement,  stress,  and  strain  by  adding  constraints, 
each  multiplied  by  Lagrange  multipliers  to  give 

1  -J  [»<*«>  -  Vi]dv-/TiV8 

v  so 

*ij  [  eij”*Js(ui, j+uj,i)J  dV  -  JY i(ui-ui*)ds  (3) 

su 

The  first  constraints,  (with  Lagrange  multiplier,  X^),  a  re  seen  to 
ensure  satisfaction  of  the  strain  displacement  equations 

eu  ■  +  “i,i>  in  v  <4) 

(where  commas  denote  partial  differentiation),  and  the  second  con¬ 
straint,  (with  multiplier,  Y^) ,  ensures 

u^  ■  Uj*  on  Su  (5) 

« 

where  Su  is  the  portion  of  the  boundary  where  displacements  are  pre¬ 
scribed,  and  u£*  are  the  prescribed  displacements  on  S u« 

In  taking  the  first  variation  of  Eq  (3) ,  not  only  are  the  dis¬ 
placements  allowed  to  vary,  but  so  too  are  the  strains  and  Lagrange 
multipliers.  It  can  easily  be  shown  (Ref  9:5)  that  the  vanishing  of 
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this  variation  requires 


where  the  n  j  are  the  direction  cosines  of  the  surface,  and  S  *  +  S^. 

If  the  are  taken  to  be  components  of  stress,  and  the  components 
of  traction,  the  first  three  volume  integrals  and  the  last  surface  in¬ 
tegral  represent  the  stress  strain  laws,  the  equations  of  equilibrium, 
the  strain  displacement  equations  and  Cauchy's  law,  respectively. 

Thus,  any  displacement  field  which  satisfies  these  equations  ensures 
that  all  but  the  last  two  surface  integrals  vanish.  The  equation  left 
to  be  satisfied  has  then  been  reduced  to 


Ji “  Ti*)6u1ds  -  ~  “  0 


These  remaining  stress  and  displacement  boundary  integrals  become  very 
important  tools  in  determining  functions  which  are  solutions  to  a  given 
problem. 
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lamic  Form  of  Reissner's  Principle 


Reissner's  principle,  as  discussed  above,  can  be  applied  only  to 
static  systems,  but  can  easily  be  extended  to  include  dynamics  with  the 
use  of  Hamilton's  principle.  To  this  end,  Eq  (3)  can  be  turned  into  a 
Lagrangian  function,  J,  by  subtracting  the  system  kinetic  energy,  K, 
represented  by 


K 


pu^u^dV 


v 


(8) 


where  p  is  the  mass  density  and  the  dots  represent  derivatives  with 
respect  to  time.  As  per  Hamilton's  principle,  (Ref  10),  the  resulting 
Lagrangian  is  integrated  over  time  from  t^  to  t£  to  give 


J 


c2 

J  { 


u^ds 


s 


a 


-  J" [  eij_1s(ui>  j+Uj  *)]  dv  “  ^yY1(ui-ui*)ds|  dt  =  0 

V  s 

u 


(9) 


This  time  dependent  version  of  Reissner's  principle  will  become  the 
basis  for  further  development  in  this  study,  and  upon  setting  the  first 
variation  to  zero,  will  be  shown  to  produce  the  equations  of  motion  and 
boundary  conditions  of  the  system  being  studied.  The  boundary  condi¬ 
tions  will  take  on  the  same  form  as  those  given  in  Eq  (7) . 


III.  Specialization  of  the  Variational  Principle 
to  Cylindrical  Shelia 

The  application  of  the  previously  described  variational  procedure 
to  a  cylindrical  shell  requires  a  slight  modification  of  Eq  (9)  to 
ensure  the  necessary  form  of  each  of  its  terms.  Before  doing  this 
however,  Eq  (9)  can  be  simplified  somewhat  by  applying  the  strain  dis¬ 
placement  equations  directly  to  the  strain  energy  expression  by  expres¬ 
sing  W  only  in  terms  of  displacements.  In  doing  so,  the  corresponding 
Lagrange  multiplier  constraint  term  on  strain  and  displacements  can  be 
eliminated  from  the  equation.  If  this  is  not  done,  an  additional 
Lagrange  multiplier  is  needed  when  Eq  (9)  is  applied  to  a  shell.  This 
is  because  for  a  shell,  W  not  only  becomes  a  function  of  mid-plane 
strains,  but  also  of  curvature.  Since  curvature  must  also  be  allowed 
virtual  changes,  this  third  set  of  Lagrange  multipliers  is  needed  to 
ensure  satisfaction  of  the  shell  curvature  displacement  relations. 
Whether  Eq  (9)  is  employed  with  W  expressed  in  terms  of  strains  and 
curvatures  in  conjunction  with  their  respective  constraint  equations,  or 
whether  W  is  put  in  terms  of  only  displacements  without  the  constraints 
is  a  matter  of  personal  preference.  In  the  preparation  of  this  study, 
both  methods  were  employed  with  identical  results.  The  latter  will  be 
described  herein. 

In  this  light,  for  a  shell  with  no  body  forces  Eq  (9)  can  be 


expressed  as 
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/!/' 


(wCu^-Js  pa^)  dv  - 


/< 


<-Mnn  64V  ul>ds 


-  f  (Xi,(3-9*)+'f1(u1-u1*)  )  ds|  dt 
Su 


(10) 


where  (i  ■  1,2,3  and  n  *  x,9) 


M  * 


nn 


3 

8* 


■  Prescribed  moments  on  S^,  positive  about  the 
positive  s  axis 

-  Prescribed  stress  resultants  on  S  defined  by: 
+  (3Mng/3s)  when  i  ■  3,  (the  z 

»  Nni  when  1^3  (the  x 

and  6  direction) 

■  Rotations  about  the  negative  s  axis 
*  Prescribed  rotations  on  S 

u 

■  Lagrange  multipliers 


direction)  and  by 


(The  subscript  i  follows  the  rules  of  the  summation  convention,  but  the 
subscript  n  does  not.)  After  taking  the  first  variation  of  Eq  (10),  the 
Lagrange  multipliers  will  form  complimentary  virtual  work  terms  similar 
to  those  in  Eq  (6).  It  follows,  the  are  tractions  and  the  Xn  is  the 
negative  of  the  moment  about  the  s  axis. 


The  application  of  Eq  (10)  to  cylindrical  shells  now  requires  the 
representation  of  the  elastic  strain  energy  density  of  the  shell  as  a 


function  of  the  displacements  of  u,  v,  w,  and  their  derivatives. 

Towards  this  end,  a  shell  coordinate  system  is  defined  as  shown  in 
Fig.  1.  For  this  system  the  thin  elastic  shell  theory  strain  displace¬ 
ment  equations  for  mid-surface  strains  reduce  to  (Ref  3) 


e  ■  u, 
x  ’x 


■  7(v,e+”) 


(ID 


Y*e  ‘ 


while  the  expressions  for  curvature  become 


K  ■  -w, 

X  XX 


(-)2(-w,fia4v,a) 


(12) 


K*e  '  I(-w'xe+i5v'x) 


The  expression  for  strain  energy  density  per  unit  area  for  an  isotropic 
cylindrical  shell  can  be  written  as  (Ref  4) 


-  *>1  [<V«e>2  - 

+  >!D2  [(K^)2  -  Id-vHK^-K^2)]  (13.) 

where 

^  -  (Eh)/ (1-v2)  (13b) 

D,  -  (Eh3)/12(l-v2) 


(13c) 


Substituting  Eqs  (11),  and  (12)  into  Eq  (13a)  produces  W  in  terms 
of  u,  v,  and  w.  Substituting  this  resultant  expression  into  Eq  (10) 
and  setting  the  first  variation  of  the  functional  to  zero  produces 


/  tU,x+av(v,0+w)  ^  6u»x+2  t(a)2(v’0+w)+aVU’x  ]5v,f 
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+  a-v)(i)2[-w,x0-tJ5v,x]9v>xjdA 
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u 
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At  this  point,  it  is  also  convenient  to  define  the  force  and  moment 
resultants  for  thin  shells  (Ref  3)  as 


Nxx  ’  Di(ex+',£e) 


"ex  ■  "xe  ■  VW-^xf, 


“ee  '  W^x5 

\z  ’  X,x^x,e> 


Mxx  ■  VV'V 

«9x  ■  «x6  ■  D2<1-',>Kx9 


«96  •  D2«8**t> 

Nez  ’  I'lVe^xe.x) 


These  are  shown  pictorially  in  Fig.  2. 

If  Eqs  (11)  and  (12)  are  substituted  into  Eqs  (15) ,  and  those 
results  placed  into  Eqs  (14) ,  the  functional  can  be  written  as 

c2 


+  IN905v’e+  i*WW 


+  Nx06v’x  +-£Nx0(Su’0)dA 
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3  3 

a  u 
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a  a  * 


To  Integrate  the  first  area  integral  in  Eq  (16),  call  it  1^,  the 
following  vectors  are  defined: 


'1  ■  N  5u  e  +  N  .<Su 
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2  "  Nx06v  ex  +  N006v  e0 


(17a) 


(17b) 


With  these 


Ix  -  yV-^jdA  +  /v.*2dA  +  /^wdA 
AAA 

-  /( N  -f^N  .  fl)5udA  -  /( N  a  +^NQa  Q)<SvdA 

J  xx,x  a  x8,0  y  x0,x  a  08,0' 


Application  of  the  divergence  theorem  gives 


■  /(|*W»  -  <8*x,x^>ce,e>s"  -  <NK0>x+iN9e,e>s'')dA 

A 

+  /(  (nxNxx4T10Nx0)6u  +  (nxNx04T19N00)6v  )ds 


where  n  are  defined  as  the  direction  cosines  between  the  normal  and 


the  y  direction.  To  integrate  the  second  integral  of  Eq  (16),  I2,  f 


scalars  are  defined  as 


vy  -  M  e  +  M  a  eQ 

X  XX  X  X0  0 


”9  ■  Mxe  ex  +  «ee 


(2 


s  6W 

X  X 


<J>Q  -  (1/ a)6w,Q 


With  Eqs  (20),  I2  can  be  written  (summing  on  y,  where  y  =*  x  or  9) 


-  /(VVVVdA 


(2 


Then,  with  the  use  of  Green's  identity. 


z2 m  "  /v3V3t i)ds 


(2 


Equation  (22)  can  be  put  in  a  more  convenient  form  by  defining 


V-<eB) 


(2 


Now,  Eq  (22)  becomes 


/*V(5w)*7x  dA  -  f  (>  Of  /3n)ds 


(2 


If  the  divergence  theorem  is  applied  to  Eq  (24),  the  result  is 


ag)6w  dA  +  /  ne(Ma6><r)Sw  ds 


6w’x(3V3n)ds 


0O*0/3n)ds 


■  y<M"e, 


cc6)6w  dA  +  ln0( MagjCt)6W  ds 


0(n«7'FQ)ds 


wliere  n  is  the  unit  normal  of  the  surface.  Finally,  after  more  inte 
gration  by  parts,  some  manipulation,  and  setting  Su^Ct^)  =  6u^(t2>  = 
Eq  (16)  can  be  written  as 


/  { y’o&.-v 


»‘Ke,x0“(a)2lSe,ee4p;')4w 
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+  y[(nxNxx"%Nxe)6u  +  (T1xNx94T10N0e)6v]  ds 
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-  yjx^g+Y^u^ds  -  y<s- B*)6xnds 


u 


-  y(ui-u1*)«Y1ds  | 


dt  =  0 


(27) 


The  three  terms  in  the  area  integral  of  Eq  (27)  can  be  recognized 
as  the  three  time-dependent  equations  of  equilibrium  for  a  cylindrical 
shell  (Ref  3:200).  The  line  integrals,  which  evolved  from  integration 
by  parts,  represent  the  boundary  conditions  which  need  to  be  satisfied 
on  the  surface  of  the  shell,  and  will  become  more  readily  recognizable 
when  applied  to  a  specific  shell  shape.  Finding  a  set  of  displacement 
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Exact  Solutions  for  a  Cylindrical  Shell 


amic  Donnell  Eauations 


Finding  a  general  solution  which  satisfies  the  equations  of  motion 
given  in  Eq  (27)  can  be  facilitated  by  following  a  procedure  used  by 
Kraus  (Ref  3:302-304).  The  first  step  is  to  apply  the  Donnell  assump¬ 
tions,  which  Donnell  proposed  in  1933  to  simplify  the  cylindrical  shell 
equations  (Ref  6) .  In  his  first  assumption,  Donnell  argued  that  the 
transverse  shearing  force  makes  a  negligible  contribution  to  the  equi¬ 
librium  of  forces  in  the  circumferential  direction  (Ref  1:200).  In  the 
second  group  of  terms  in  the  area  integral  of  Eq  (27) ,  the  terms  which 
make  up  this  force  are 


N  ■  — fM  +  )  at  0 

8z  a'  x0,x  +  a  08,0  u 


and  can  thus  be  eliminated.  Substitution  of  Eqs  (11),  (12)  and  (15) 
into  the  equations  of  motion  from  Eq  (27)  (given  in  the  first  area 
integral  after  this  first  Donnell  assumption  has  been  made) ,  allows  the 
equations  to  be  represented  in  terms  of  u,  v,  and  w.  Through  a  series 
of  operations  and  substitutions,  the  following  three  new  equations  can 
be  derived 


v\u)  +  -  (i)3w,x6e 
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xxxx 


2(1-Ki)  32 


1_  f /l-v2  32  3-v?2\ 

t2  Lr**  at2 "  / 


^^pw  +  (|)2W  -I-  (h2/12)74(w)^ 


+  >s(l-u)V4(w)  +  (Vw.  +  ( 

a  X-X 


(i)4w'ee] 


(28b) 


(28c) 


These  new  equations,  which  are  commonly  known  as  the  Donnell  equations, 
are  of  higher  order  than  the  original  system  of  equations,  (one  eighth 
order  plus  two  fourth  order  as  compared  to  one  fourth  order  plus  two 
second  order) .  The  advantage  to  this  increase  in  order  is  that  the  dis¬ 
placements  have  become  decoupled.  However,  any  solution  that  is  found 
to  satisfy  these  new  equations  should  also  be  checked  back  in  the  ori¬ 
ginal  system  to  guard  against  the  introduction  of  extraneous  roots. 

Solution  Method  for  Closed  Cylindrical  Shells 

Kraus  (Ref  3)  outlines  a  method  for  determining  the  natural  fre¬ 
quencies  for  a  closed  cylindrical  shell  by  first  assuming  a  solution 
function  given  by  (symmetric  modes) 
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V  ” 


£ 


Biexp(Xix/l)8in(n9)cos  (tut) 


w 


K 

53  Ciexp(Xix/l)cos(n6)cos(<jdt) 


(29) 


where  A^,  B^,  and  are  arbitrary  constants,  and  the  X^  are  to  be 
determined  so  that  Eqs  (29)  satisfy  the  governing  equations.  Due  to  the 
order  of  the  system,  the  number  of  terms  in  the  summation,  K,  can  be 
shown  to  be  eight  (Ref  7) .  If  these  assumed  solutions  are  substituted 
into  the  dynamic  Donnell  equations,  Eqs  (28),  the  result  is 


Aif2^ 

CiLi-v 


(3-v)An2y 

Pv 


+ 


+  n2  (I+vZ)l 


(30) 


*  2a3  -  A2  |^2+(3-\))n2Y  +  (2/5)n4Y2J 
+  A|(3-v)n2Y  -  2n2(l-v2Z)  +  (l-v)n4y2 

+  (l/C)(3-v)n6Y3]  -  (l/£)(l-v)n8Y4  (32) 

where 

A  -  pw2a2(l-v2)/E  (33) 

C  -  12(a/h)2  (34) 

Z  -  (Xia)2/(nl)2 

and 

Y  =*  (1-Z)  (35) 

At  this  point,  the  solution  to  a  closed  cylindrical  shell  becomes  fairly 
straight-forward.  There  are  four  boundary  conditions  to  be  satisfied  at 
each  end.  They  are  (Ref  1:136) 


N 

XX 

-  0 

or 

u  ■  0 

(36) 

M 

XX 

-  0 

or 

"■*  ■ 0 

(37) 

-  0 

or 

v  ■  0 

(38) 

[WT(Mne,9) 

-  0 

or 

w  ■  0 

(39) 

If  Eqs  (29)  are  substituted  into  the  equations  for  the  eight  boundary 
conditions,  the  result  is  eight  homogeneous  equations  in  the  eight  un¬ 
knowns  C^,  (A^  and  can  be  eliminated  with  Eqs  (30)  and  (31))  •  The 


solution  to  these  eight  equations  is  found  by  setting  the  determinant 
of  their  coefficients  to  zero.  The  determinant  is  dependent  upon  the 
X^,  which  are  functions  of  uj  and  n  through  the  eighth  order  Eq  (32). 
Thus,  for  a  given  n,  and  ui  is  chosen  and  Eq  (32)  solved  for  the  eight 
X^.  The  determinant  is  then  calculated  with  these  values.  The  pro¬ 
cedure  is  repeated  in  an  iterative  manner  until  an  u  is  found  where  the 
determinant  vanishes. 

As  a  point  of  interest,  examination  of  Eq  (32)  shows  that  the 
eighth  order  equation  in  X^  has  only  even  powers.  This  leads  to  a  dis¬ 
tinct  relation  between  the  eight  roots;  four  roots  are  the  negatives  of 
the  remaining  four. 

Unfortunately,  the  above  method  cannot  be  used  for  open  cylindrical 
shells  because  four  more  boundary  conditions  are  introduced  on  each  of 
the  edges  of  constant  0.  They  are  (Ref  1:157) 


Nxe 

-  0 

or 

u  ■  0 

(40) 

Nee 

-  0 

or 

v  -  0 

(41) 

^z^xQ.x 

-  0 

or 

w  ■  0 

(42) 

Mee 

-  0 

or 

« 

CD 

a 

o 

(43) 

There  are  now  a  total  of  sixteen  boundary  conditions  to  be  satisfied 
and  only  eight  undetermined  coefficients.  Thus,  Eqs  (29)  are  not  gene¬ 
ral  enough  to  provide  a  solution  for  any  one  n.  For  this  reason,  the 
variational  procedure  described  earlier  will  be  used  in  an  attempt  to 
find  an  approximate  solution  by  superposition  of  functions  with  several 


values  of  n. 


As  discussed  in  the  previous  section,  the  presence  of  eight  addi¬ 
tional  boundary  conditions  for  the  open  cylindrical  shell  prohibits  the 
use  of  the  procedure  outlined  by  Kraus,  (Ref  3),  to  find  a  solution  for 


the  free  vibration  modes.  To  overcome  this  obstacle,  the  variational 
procedure  developed  earlier  and  expressed  by  Eq  (27)  will  be  used.  The 
approach  will  be  to  find  a  trial  solution  function  which  satisfies  Eq 
(27) ,  and  hence  satisfies  or  approximates  the  governing  shell  equations 
and  all  sixteen  boundary  conditions. 

Determination  of  Surface  Integrals 

The  boundary  Integrals  of  Eq  (27)  can  now  be  expressed  in  terms  of 
a  specific  shell  shape.  Thus  far,  Eq  (27)  applies  to  all  cylindrical 
shells  with  the  only  assumption  being  the  thickness,  h,  is  small  as  com¬ 
pared  to  the  radius,  r.  Specializing  Eq  (27)  to  an  open  cantilevered 
shell,  such  as  the  one  shown  in  Fig.  3,  allows  the  surface  integrals  to 
be  written  explicitly.  Before  doing  so  however,  the  prescribed  force 
and  moment  resultants  on  Sff,  and  the  prescribed  displacements  on  are 
all  set  to  zero.  This  eliminates  the  terms  with  the  asterisk  in  the 
last  four  integrals  of  Eq  (27) .  The  remaining  integrals  from  Eq  (27) , 
if  written  such  that  the  surface  integrations  will  be  done  counter¬ 
clockwise  (as  viewed  from  the  +z  direction)  around  the  four  edges  with 
the  sign  convention  of  Fig.  2,  become 
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(44) 


Several  terms  in  this  equation  can  now  be  eliminated  by  inducing 
the  second  Donnell  assumption  (Ref  3) .  It  can  be  reasoned  the  changes 


in  the  stretching  displacement  in  the  circumferential  direction,  v, 
have  little  effect  on  the  curvature  and  twist  of  the  mid-section. 

Looking  at  Eqs  (12) ,  this  assumption  is  tantamount  to  saying  derivatives 
of  v  are  small  compared  to  derivatives  of  w.  This  allows  the  v/a  terms 
to  be  eliminated  in  the  eighth  boundary  integral  in  Eq  (44) .  With  both 
Donnell  assumptions,  the  area  integral  in  Eq  (44)  represents  a  new  set 
of  equations  which  again  can  be  transformed  into  the  Donnell  equations, 
Eqs  (28) ,  and  the  boundary  terms  become  the  more  recognizable  shell 
boundary  conditions  for  an  open  cylindrical  shell  given  in  Eqs  (36) 
through  (43) . 


Determining  a  Trial  Solution  Function 


In  selecting  a  trial  solution  function  which  is  to  be  used  to 
satisfy  Eq  (44)  and  hence  be  a  solution  to  the  open  shell,  it  is  advan¬ 
tageous  to  start  with  one  which  satisfies  as  much  of  Eq  (44)  term  by 
term  as  possible.  In  doing  so,  those  parts  of  Eq  (44)  which  are  satis¬ 
fied  in  this  manner  can  be  eliminated  and  the  complexity  of  the  calcu¬ 
lations  reduced.  Therefore,  if  Eqs  (29)  are  used  to  build  a  trial 
function,  then  the  area  integral  in  Eq  (44)  may  be  dropped.  This  is 
because,  as  long  as  Eqs  (30),  (31)  and  (32)  are  enforced,  Eqs  (29)  do  in 
fact  satisfy  the  differential  equations  of  motion. 

The  flexibility  of  the  variational  procedure  allows  a  trial  solu¬ 
tion  to  be  used  which  does  not  have  to  meet  the  required  boundary  condi¬ 
tions  term  by  term.  This  feature  will  be  employed  by  using  Eqs  (29) . 

For  example,  if  Eqs  (29)  are  used  to  satisfy  the  necessary  boundary  con¬ 
dition,  Nqq  -  0,  on  an  edge  of  constant  0,  then  the  equation  that  must 
be  satisfied  is 

fexpdjx/l)  (B^n/a+C^/a+vAjXj/1)  ]  cos(±<=n)  *  NQ0  ~  0  (45) 

For  N00  to  be  zero  for  any  x  along  a  free  edge,  the  term  outside  the 
brackets  must  vanish.  This  requires  that  n*(2m+l)Tr/2«.  If  this  result 
is  placed  into  Eqs  (29) ,  the  displacements  u  and  w  are  found  to  vanish 
along  the  free  edges  of  constant  0.  Since  this  is  totally  undesirable 
for  the  cantilevered  shell,  the  variational  procedure  must  ensure  satis¬ 
faction  of  the  boundary  conditions  on  these  edges  through  Eq  (44). 

At  this  point,  Eqs  (30),  (31),  and  (32)  can  be  enforced  on  Eqs  (29), 
and  the  result  written 
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u  -  ^  4'^C^exp(X^x/l)cos(n8)cos(u)t) 

i=*l 

8 

v  -  £  $.C.exp(X.x/l)sin(n0)cos(uit)  (46) 

i-1 

8 

w  »  £  C.exp(X.x/l)cos(n6)cos(ajt) 

i-1 

where  ^  and  ^  are  determined  from  Eq  (30)  and  (31)  utilizing  the  X^ 
determined  from  Eq  (32). 

Thus  far,  Eqs  (46)  can  be  considered  each  one  "term"  consisting  of 
eight  parts  through  the  constants  summed  on  i  from  one  to  eight.  If 
seven  boundary  conditions  are  chosen  to  be  satisfied  by  Eqs  (46)  ,  seven 
of  the  C's  can  be  written  in  terms  of  the  eighth.  Performing  this 
operation  also  tailors  the  trial  solution  functions  to  be  formed  as  a 
series  of  N  terms.  Each  term  is  of  the  form  of  Eqs  (46),  so  it  satis¬ 
fies  the  seven  chosen  boundary  conditions  exactly,  term  by  term.  Hence, 
those  corresponding  integrals  in  Eq  (44)  may  be  eliminated.  The  seven 
boundary  conditions  that  are  to  be  used  must  be  chosen  from  Eqs  (36) 
through  (39).  Equations  (40)  through  (43)  cannot  be  used,  because  as 
discussed  above,  Eqs  (46)  cannot  be  used  to  satisfy  boundary  conditions 
both  along  the  edges  of  constant  0  and  of  constant  x.  A  choice  must  now 
be  made  as  to  which  of  these  seven  boundary  conditions  will  be  satisfied 
exactly,  and  which  one  is  to  be  left  to  be  approximated  through  the 
variational  technique.  For  some  problems  there  are  boundary  conditions 
that  can  be  considered  to  be  more  crucial  than  others.  For  a 


cantilevered  shell,  not  requiring  that  the  w  displacement  at  the  fixed 
end  be  zero  would  lead  to  erroneous  results  in  the  natural  frequencies 
and  mode  shapes.  However,  a  v  displacement  which  is  left  arbitrary  at 
the  fixed  end  represents  a  shell  supported  in  a  channel  which  restricts 
motion  in  all  but  the  surface  tangential  direction  (the  v  direction) . 
Such  a  shell  could  be  expected  to  have  natural  frequencies  and  mode 
shapes  very  similar  to  those  of  a  shell  with  the  end  totally  clamped. 
Therefore,  the  boundary  condition  v-0  at  x-0  will  be  left  to  be  approxi¬ 
mately  satisfied  by  the  variational  principle,  rather  than  exactly. 

Equations  (46)  are  now  substituted  into  the  necessary  seven 
boundary  conditions  in  Eqs  (36)  through  (39)  to  set  the  three  remaining 
displacements  to  zero  at  x-0  and  the  four  force  and  moment  resultants  to 
zero  at  x=l.  This  produces  seven  equations  in  the  eight  constants  C^. 
After  some  manipulation  the  seven  equations  can  be  written 

If  u(o,0)  -  0,  then 

If  w(o,9)  *  0,  then 

8 

If  w,x(o,9)  -  0,  then  £  XiCi  *  0  (47c) 

i-1 

8 

If  M^U.9)  -  0,  then  12  exp(X1)C1  [-(A^l)2  +  v(n/a)2J- 


8 

£  **C.  =  0  (47a) 

i-1  1 


8 

£  Ct  -  o  (47b) 

i-1 


0  (47d) 


If  N  (1,0)  -  0,  then  53exP(Xi)Ci|4'iXi/1  +  Iv(n<f,i+1)J  “  0  (47e) 

505  i-1  L 

If  Q  (1,0)  -0,  then  23exp(Xi)C.Xi f(Xi/l)2-(2-v) (n/a)2J  -  0  (47f) 

xz  i-I  L 

If  **  -  o,  then  Eexp(xi)ci[Di(4>ixi/;L  ■  iTl4'i) 

a  ^  ^ 

+  (^)2D2(l-v)Xin/l]  -  0  (47g) 

Using  these  seven  equations  and  solving  for  C2  through  Cg  in  terms  of 
gives 

c±  -  ric1  (47h> 

Since  the  \±  are  complex,  and  r±  are  also  complex.  For  this 

reason,  Eq  (47h)  cannot  be  solved  by  using  a  standard  linear  equation 
solution  technique  such  a  Gaussian  reduction.  Finding  the  inverse  of 
the  solution  matrix  is  one  possible  way  of  finding  I^. 

With  this  expression  Eqs  (46)  can  be  written 


u  -  ci  [ri'i'iexp(Xix/l)cos(n0)cos(aJt)] 

v  -  C1[ri$iexp(X1x/l)sin(n0)cos(o,t)]  (48) 

w  -  C1[riexp(X1x/l)cos(n0)cos(u>t)] 

The  summation  on  i  from  this  point  on  will  be  assumed  to  be  from  one  to 


eight,  and  the  summation  sign  dropped.  Therefore,  these  three  dis¬ 
placement  functions  will  be  considered  to  be  one  "term"  consisting  of 
eight  parts.  The  only  unknown  constant  in  Eqs  (48)  is  C^.  Eqs  (21) 
are  solutions  to  the  three  shell  differential  equations,  and  all  but 
one,  (v«0  at  x*0),  of  the  boundary  conditions  at  the  edges  of  constant 
x.  The  constants  and  X^  are  determined  for  a  given  a>  and  n 

from  Eqs  (30),  (31),  and  (32).  The  components  of  are  determined 
from  the  seven  boundary  conditions  as  described  above. 

The  assumed  solution  functions  needed  to  satisfy  Eq  (44)  can  now 
be  defined  by  creating  a  series  using  Eqs  (48) .  Each  function  consists 
of  N  terms  given  by 

N 

u  *  2  E  r.4'.exp(X.x/l)cos(n0)cos(u)t) 
n-0 

N 

v  ”  £  En  ri<|»iexp(X1x/l)cos(n0)cos((Dt)  (49) 

n**0 

N 

w  -  £e  r .exp ( X ±x/ 1 ) cos (n0 ) cos (u> t) 
n-0 

Each  "term"  in  these  series  consists  of  eight  separate  parts  summed  on 
i.  The  same  constants,  E^,  appear  in  each  of  the  three  functions,  as 
did  C^,  because  the  differences  between  the  ratios  of  u,  v,  and  w  are 
accounted  for  in  the  and  Thus,  Eqs  (49)  provide  a  solution 

function  where  each  term  satisfies  exactly  the  area  integral  and  seven 
of  the  boundary  integrals  of  Eq  (44).  Also,  the  parameter  n,  which 
could  have  had  any  value  up  to  this  point,  has  been  restricted  to 


Integer  values.  This  choice  Is  acceptable  as  long  as  n  is  chosen  so 
as  to  Include  an  acceptable  class  of  functions  for  utilization  in  the 
variational  procedure.  In  the  assumed  solution  form,  Eq  (29),  n  repre¬ 
sents  the  number  of  longitudinal  nodal  lines  in  the  mode  shape  for  a 
given  natural  frequency.  Thus,  n  should  take  on  values  of  0  through  N 
to  provide  a  "complete"  set  of  functions. 


Obtaining  a  Solution 

Equation  (44)  used  in  conjunction  with  Eqs  (49)  can  now  be  used  to 
solve  for  the  desired  vibrational  modes  for  the  open  cantilevered 
cylindrical  shell.  As  discussed  above,  the  terms  which  are  satisfied 
identically  by  Eqs  (49)  can  be  eliminated  from  Eq  (44)  .  This  leaves 
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as  the  equation  needed  to  be  satisfied  to  ensure  the  functional  is  at  a 


stationary  value 


The  general  form  of  each  term  in  Eq  (50)  is  a  force  multiplying  a 
virtual  displacement,  or  a  displacement  multiplying  a  virtual  force. 

In  this  light,  each  term  represents  a  form  of  virtual  work.  The  func¬ 
tions  in  Eqs  (49)  are  complex,  but  the  real  and  imaginary  parts  each 
Independently  satisfy  the  differential  equations.  Only  the  real  parts 
should  be  used  when  they  are  substituted  into  Eq  (50) .  In  this  manner 
each  term  is  a  real  force  multiplying  a  real  displacement  and  takes  the 
form 
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I  If' 


Re(F)  Re(5u) 


dt 


resulting  in  a  real  work  term.  Here  F  is  the  generalized  complex  force 
and  u  is  the  generalized  complex  displacement.  If  the  integration  with 
respect  to  time  is  taken  to  be  one  period,  it  can  easily  be  shown  that 
the  same  results  can  be  achieved  by 


J  *sRe(F6u)  ds 
s 

where  u  is  the  complex  conjugate  of  u.  This  method  will  be  used  when 

Eqs  (49)  are  substituted  into  Eq  (50) .  Also  in  performing  the  substitu 

tion  the  most  general  arbitrary  variation  of  a  displacement  or  force 

resultant  is  taken  to  be  a  variation  on  the  coefficients  E  .  For  ex- 

n 

ample,  substitution  of  v  from  Eq  (49)  into  the  5v  in  the  first  term  of 
Eq  (50)  would  result  in 
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After  substituting  in  Eqs  (49)  in  this  manner  and  combining  terms,  Eq 


(50)  can  be  written 


sin(n»)cos(m*)  ^riexp(Xix/l)  (<J>iXi/l-A'i)J 
[^exp^x/l)  J  5Emdx 

1 

+  J'lE^p^  cos(n«)sin(ni“)  ^riexp(X^x/l)  (-^^+^4v4'iXi/l)J 
0 

^ jrjexP(XjX/i) J  5Emdx 
l 

+  f ^2vE sin(n*)cos(m*)  ^I^expCX^x/l) (2-v) (X^/1)2-(^)2J 
jexp(XjX/l) J  SE^dx 

1 

+ J'^2E^p  cos(n*)sin(m“)  ^r^exp^jx/l)  ((^-)2-\>(Xi/l)2)J 
0 

^mrjexp(XjX/l)J  <5Emdx 

cc 

+  /En  [ri*l]  [Vj^jV1  -  ?”V 

— oc 

+  (— )2D(l-v)T .mX ./I  I  sin(n0)sin(m6)  ad0 


*  0  (52) 

Since  this  expression  must  be  valid  for  arbitrary  values  of  <5Em>  it  is 
convenient  to  make  the  choice 
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(53) 
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With  this  choice  Eq  (52)  becomes  a  system  of  M+l  equations  each  with 
N+l  terms  multiplying  the  constants  Em«  Setting  N*M,  Equation  (52)  can 
be  represented  by  the  complex  matrix  equation 

E^  -  0  n,m  -  0,1,2...,N  (54) 

Each  term  in  the  square  matrix  is  found  by  solving  Eq  (52)  for  the 
appropriate  values  of  m  and  n  for  each  coefficient  to  E^.  It  is  very 
Important  to  remember  each  Integral  consists  of  a  series  of  eight  terms 
summed  on  1  multiplying  a  series  of  eight  terms  summed  on  j .  Even 
though  this  integration  of  the  product  of  two  series  seems  laborious, 
the  integrals  of  Eq  (52)  reduce  to  simple  summations  which  are  handled 
very  nicely  on  a  digital  computer.  The  solution  to  the  set  of  equations 
represented  in  Eq  (54)  can  be  found  by  setting  the  determinant  of  Q  to 
zero.  The  "eigenvalue"  of  the  system  is  w  and  is  implemented  through 
Eq  (32).  Repeatedly  choosing  values  for  u>,  performing  the  necessary 
calculations  to  fill  Q,  and  then  computing  its  complex  determinant  pro¬ 
vides  a  procedure  which  may  be  used  to  find  zeroes  of  det(Q)  iteratively. 
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and  thus  solve  for  approximate  fundamental  frequencies  of  the  shell. 

In  general,  when  working  with  structures  which  have  linear  force- 
displacement  relationships,  it  can  be  shown  through  the  Betti-Maxwell 
reciprocal  theorem  that  the  system  matrix  is  symmetric.  It  is  expected 
therefore,  that  this  system  derived  through  the  variational  procedure 
should  not  be  an  exception  (See  Appendix  A).  For  complex  matrices,  the 
analog  to  symmetric  is  the  Hermitian  form.  It  is  not  readily  apparent 
from  studying  Eq  (52)  that  Q  is  Hermitian.  This  topic  will  be  discussed 
more  in  the  results  section. 


VI.  Analysis  of  Results 


The  natural  frequencies  and  mode  shapes  of  a  cantilevered  open 
cylindrical  shell  were  determined  experimentally  with  a  Hewlett  Packard 
5423  Structural  Dynamics  Analyzer.  The  dynamics  analyzer  can  accurately 
determine  the  natural  frequencies  and  mode  shapes  of  structures  such  as 
a  cantilevered  shell,  and  thus  provide  an  excellent  verification  of 
analytic  results. 

Experimental  Procedures 

The  HP  Structural  Dynamics  Analyzer  is  a  two  channel  fast  Fourier 
series  spectrum  analyzer  specifically  designed  for  the  purpose  of  cal¬ 
culating  mode  shapes  and  natural  frequencies  of  structures.  The  analy¬ 
zer  performs  this  function  utilizing  two  input  signals  generated  from 
impact  testing.  The  signals  are  sent  to  the  analyzer  in  the  form  of  (1) 
a  structure  input  signal  from  a  force  transducer  providing  the  impact 
time  history  and  (2)  a  structural  output,  or  response  time  history,  from 
an  accelerometer  measuring  the  structure  response  to  the  impact.  A 
sample  plot  for  each  of  these  signals  is  given  in  Fig.  4.  This  figure 
is  representative  of  the  response  at  one  point  due  to  an  impact  at  ano¬ 
ther.  The  input  signal  time  history  is  clearly  depicted  as  a  force  im¬ 
pulse,  while  the  output  shows  the  resultant  vibration  response.  A 
transformed  ratio  of  this  output  signal  to  the  input  signal  provides  a 
transfer  function  for  the  respective  two  points.  Figure  5  represents 
the  transfer  function  of  the  measurements  given  in  Fig.  4.  The  spikes, 
or  commonly  called  poles,  in  the  plot  at  various  frequencies  are  a  re¬ 
sult  of  the  transfer  function’s  denominator  approaching  zero,  causing  a 


large  Increase  In  the  value  of  the  function  itself.  Synonymous  with 
the  Laplace  transform  representation  of  the  transfer  function,  these 
poles  occur  at  the  natural  frequencies  of  the  system  being  analyzed. 
Thus,  by  locating  these  poles  the  natural  frequencies  can  be  calculated. 
This  transfer  function  provides  one  element  of  the  transfer  function 
matrix  G.  Each  element,  for  example,  represents  the  ratio  of  the 
output  at  point  j  due  to  the  input  at  point  i.  If  a  structure  is  de¬ 
fined  by  N  points,  a  complete  transfer  function  matrix  is  defined  as  an 
N  by  N  matrix.  If  enough  impact  data  is  taken  to  fill  one  row  or  one 
column  of  this  matrix,  the  dynamics  analyzer  can  generate  the  mode 
shapes  and  natural  frequencies  of  the  structure.  (One  row  or  column 
provides  a  minimum  amount  of  data,  but  more  can  be  generated  if  a 
greater  degree  of  accuracy  is  desired.) 

The  experiment  was  performed  utilizing  the  shallow  cantilevered 
shell  depicted  in  Figures  6a  and  6b.  The  shell  as  shown  in  Fig.  6b  was 
defined  to  the  analyzer  to  consist  of  56  points.  These  points  were  cho¬ 
sen  to  provide  a  "fine"  enough  representation  of  the  structure,  while  at 
the  same  time  keeping  the  number  of  points  at  which  data  must  be  taken 
to  a  minimum.  The  shell  dimensions  were  chosen  because  some  limited 
numerical  results  were  found  for  this  specific  shape  (Ref  5) .  These 
numerical  results  were  generated  using  the  Ritz  method  and  provide  as¬ 
surance  that  the  experimental  setup  and  methods  are  accurate. 

The  shell  was  rigidly  clamped  on  a  fixed  support  resulting  in 
points  8,  9,  24,  25,  40,  41  and  56  being  fixed.  After  an  impact  point 
was  chosen,  the  output  response  was  measured  in  the  radial  direction  at 
all  56  points,  and  in  the  tangential  direction  at  all  but  the  centerline 


Figure  6b.  Top  View  of  Modal  Analyzer  Shell 


and  fixed  points.  The  centerline  points  were  omitted  because  erroneous 
motion  could  be  generated  in  the  mode  shapes  if  the  accelerometer  was 
not  placed  exactly  on  the  centerline.  Having  both  of  these  measure¬ 
ments  allows  the  w  and  v  displacements  to  be  included  in  the  mode  shapes, 
thus  making  them  more  accurate,  but  they  do  not  affect  the  frequency 
results.  As  a  side  note,  the  accuracy  of  the  results  was  dependent  upon 
how  well  the  impact  point  was  chosen.  If  the  impact  point  happened  to 
lie  on  a  nodal  line  of  one  of  the  normal  modes,  then  that  mode  was  not 
excited  sufficiently.  Therefore,  several  trial  runs  must  be  performed, 
or  the  characteristics  of  the  mode  shapes  must  be  known  well  enough  a 
priori  in  order  to  find  an  acceptable  excitation  point.  For  this  test, 
point  6  or  54  was  found  to  work  very  well. 

Experimental  Results 

Table  I  compares  the  experimental  results  to  those  calculated  using 
the  Ritz  method  (Ref  5).  The  results  can  be  seen  to  compare  favorably. 
The  last  column  of  the  table  shows  that  the  experimental  data  is  between 
five  and  ten  percent  lower  than  the  analytic  results.  The  experimental 
data  is  expected  to  be  lower  due  to  the  difficulty  in  achieving  a  true 
clamped  boundary  condition,  and  the  additional  mass  of  the  accelerome¬ 
ter.  The  mode  shapes  for  these  natural  frequencies  (Ref  5)  are  presen¬ 
ted  in  Fig.  7.  From  the  figure  and  Table  I  it  can  be  seen  that  the 
first  mode  of  vibration  for  the  shell  is  a  twisting  mode,  while  the 
second  is  a  bending  mode.  This  is  unlike  the  vibration  characteristics 
of  a  cantilevered  plate,  which  has  a  symmetric  bending  mode  for  its 
first  fundamental  frequency,  and  a  twisting  mode  for  its  second  (Ref 
12).  Figures  8  and  9  show  the  shell's  normalized  mode  shapes  as 
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calculated  by  the  modal  analyzer.  Figure  8  clearly  shows  the  antisym¬ 
metric  twisting  motion,  while  Fig.  9  shows  an  antisymmetric  motion  in 
the  circumferential  direction.  Figure  8  also  reveals  a  reverse  bending 
of  the  tips  of  the  shell  at  the  two  free  corners.  Since  these  two  areas 
undergo  greater  amplitude  and  hence  experience  greater  accelerations 
than  the  other  regions  of  the  shell,  they  are  more  susceptible  to  the 
increased  mass  of  the  attached  accelerometer.  This  increase  in  mass 
would  cause  the  corners  to  "lag"  behind  the  normal  motion  of  the  shell 
as  shown  in  the  figure.  Also,  in  Fig.  9  a  noticeable  difference  exists 
between  the  motion  on  the  left  edge  of  the  shell  versus  the  right  edge. 
The  left  edge  happens  to  be  the  side  where  the  impact  point  was  chosen 
(point  6) ,  and  therefore  should  record  a  slightly  greater  response  to 
the  Impact  than  the  far  side.  This,  along  with  the  fact  that  the  motion 
in  Fig.  9  has  been  magnified  several  times  to  make  it  visible,  help  to 
explain  this  phenomenon,  which  was  also  present  in  the  other  modes.  (See 
Appendix  B  for  the  remaining  mode  shapes.) 

All  of  the  mode  shapes  calculated  by  the  modal  analyzer  compared 
identically  to  those  calculated  with  the  Rltz  technique  with  the  excep¬ 
tion  of  the  seventh  mode.  This  mode  shape,  given  in  Figures  10  and  11, 
does  not  agree  with  the  one  given  in  Fig.  7.  An  investigation  into  this 
anomaly  revealed  that  the  seventh  mode,  as  determined  by  the  analyzer, 
is  actually  a  combination  of  the  seventh  and  eighth  modes.  Table  1 

I 

shows  that  the  natural  frequencies  of  these  two  modes  are  very  close 
together.  In  addition.  Fig.  7  shows  that  the  two  mode  shapes  have  very 
similar  nodal  lines,  but  that  mode  seven  is  antisymmetric  while  eight  is 
symmetric.  If  the  two  mode  shapes  are  superimposed,  the  displacements 


on  one  side  would  cancel  each  other  while  on  the  opposite  side  they 
would  be  additive.  This  is  exactly  what  is  depicted  in  Fig.  10.  Thus, 
to  separate  these  modes  on  the  analyzer,  each  one  must  be  individually 
excited,  while  the  other  is  suppressed.  This  could  be  done  by  very 
carefully  choosing  the  Impact  point,  or  more  easily,  by  exciting  the 
structure  with  a  shaker  at  a  specific  frequency  rather  than  by  an  impact 
(see  Ref  13) . 

Computational  Method  Development 

The  variational  procedure  derived  earlier  and  represented  by  Eq 
(52) ,  was  used  as  a  basis  for  the  development  of  a  computer  program  to 
be  used  to  calculate  the  vibration  characteristics  of  a  cantilevered 
shell.  The  computer  program  was  written  to  perform  the  necessary  opera¬ 
tions  to  fill  the  Q  matrix  and  solve  for  its  determinant.  The  roots  of 
the  equation,  and  hence  the  shell’s  natural  frequencies  are  found  in  an 
iterative  manner  by  finding  values  of  u>  which  cause  the  determinant  to 
vanish.  The  computer  program  consists  of  a  straight-forward  coding  of 
Eqs  (30),  (31),  (32)  and  (47h)  to  solve  for  the  components  4^, 
and  These  are  then  used  in  conjunction  with  Eq  (52)  to  fill  the  Q 

matrix.  All  of  the  required  operations  were  easily  performed  by  inter¬ 
nal  or  system  subroutines.  A  brief  description  and  a  listing  of  the 
program  is  given  in  Appendix  C.  The  program,  written  in  Fortran  V,  can 
fill  and  calculate  the  determinant  of  a  ten  term  matrix  in  approximately 
one  second  of  computation  time  on  a  CDC  6600  digital  computer  (NOS  oper¬ 
ating  system) . 

The  computer  program  was  developed  up  to  the  point  of  performing 
all  of  the  necessary  calculations  to  solve  for  the  natural  frequencies. 


However,  the  initial  program  results  indicate  that  further  work  must  be 
done  to  improve  the  accuracy  of  the  program's  calculations.  Simple 
tests  show  that  numerical  problems  are  present  and  are  significantly  af¬ 
fecting  the  results.  For  instance,  after  solving  Eq  (32)  for  its  eight 
roots,  one  of  them  must  be  established  as  X^.  While  this  choice  is  com¬ 
pletely  arbitrary  at  first,  the  same  root  must  then  be  chosen  each  suc¬ 
cessive  iteration  to  ensure  that  the  determinant  of  Q  behaves  as  a  con¬ 
tinuous  function.  While  the  choice  significantly  affects  the  magnitudes 
of  the  r.,  and  hence  the  elements  of  Q  itself,  it  should  not  change  the 
solutions  to  Eq  (52).  Unfortunately,  this  indeed  occurs,  and  different 
solutions  can  be  found  depending  upon  how  is  chosen.  A  way  to  elimi¬ 
nate  this  problem  must  be  found  before  accurate  solutions  to  the  problem 
can  be  expected. 

As  mentioned  earlier,  the  matrix  Q  in  Eq  (54)  es  expected  to  be 
Hermitlan,  which  would  result  in  the  determinant  of  Q  being  purely  real. 
Investigating  the  elements  of  Q  for  different  values  of  oa  shows  that  Q 
is  not  Hermitlan,  and  that  its  determinant  is  complex.  However,  the 
matrix  diagonal  is  purely  real  (another  requirement  of  a  Hermitlan 
matrix) ,  and  some  specific  elements  are  very  close  to  being  the  complex 
conjugates  of  one  other.  Again,  the  numerical  problems  discussed  above 
must  be  eliminated  before  any  further  analysis  can  be  performed  in  this 


area 


VII.  Conclusions  and  Recommendations 

The  variational  procedure  which  has  already  been  used  to  calculate 
the  vibrational  characteristics  of  a  membrane  and  plate  with  mixed 
boundary  conditions  can  be  extended  to  a  cylindrical  shell.  A  computa¬ 
tional  method  has  been  derived  to  perform  all  of  the  necessary  opera¬ 
tions  needed  to  implement  the  procedure,  but  the  nature  of  the  calcula¬ 
tions  makes  them  sensitive  to  numerical  errors. 

The  type  of  problem  which  may  be  solved  with  the  variational  tech¬ 
nique  can  easily  be  extended  from  the  cantilevered  open  cylindrical 
shell  that  was  presented  in  this  study.  Any  form  of  mixed  boundary  con¬ 
ditions  can  be  implemented  by  changing  the  boundary  integrals  that  were 
evaluated  around  the  four  edges.  In  this  manner,  virtually  any  combina¬ 
tion  of  boundary  conditions  can  be  investigated.  The  technique  can  even 
be  modified  to  analyze  shells  with  edge  cracks.  For  this  case,  the 
boundary  Integrals  must  be  evaluated  along  two  additional  free  edges. 

In  this  light,  the  variational  technique  is  a  powerful  method  which 
can  be  applied  to  problems  which  cannot  be  handled  by  techniques  such  as 
the  Rltz  method  and  the  Galerkin  method.  Since  this  flexibility  would 
be  applicable  to  a  wide  range  of  practical  problems,  further  work  in 
improving  the  accuracy  of  the  numerical  computations  is  justified. 
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13. 


The  BetCi  and  Rayleigh  Reciprocal  Theorem  states  that  if  a  body  is 
subjected  to  two  systems  of  body  and  surface  forces,  then  the  work  that 


would  be  done  by  the  first  system  T^,  in  acting  through  the  displace 
ments  due  to  the  second  system  of  forces  is  equal  to  the  work  that 
would  be  done  by  the  second  system  T^',  in  acting  through  the  dis¬ 
placements  u^  due  to  the  first  system  of  forces  (Ref  11:391).  For  a 
system  with  no  body  forces  this  can  be  expressed  as 


J * T^u^"  ds  “  J T^'  u^  ds 
s  s 

if  solutions  are  assumed  to  be  of  the  form 

T±  -  ReCl^expdut)) 

^  -  Re(unexp(iu>t)) 


T±'  -  ReCT^expUwt)) 


u^'  -  Re(umexp(iut)) 


and  are  substituted  into  Eq  (A-l),  the  result  is 

J“  ReC^expCiut))  Re(umexp(iu>t) )  ds 

s 

*  ^yReCl^expdwt) )  Re(unexp(io»t))  ds 


(A-l) 


(A-2) 


(A-3) 


This  can  be  written  as 


^*[Re(Tn)cosa)t  -  Im(Tn)sinuit]  [Re(um)cosa)t  +  Im(um)siniot]  ds 

3 

-  [Re(Tm)cos<»)t  -  Im(Tm)sin£i)t]  [Re(un) cosoit  -  Im(un)sinu)t]  ds  (A-4) 


If  Hamilton's  principle  is  applied,  and  these  two  integrals  are  inte¬ 
grated  with  respect  to  time  over  one  period,  the  result  is 


h  J"  [Re(Tn)Re(um)  +  Im(Tn)Im(um) ]  ds 
s 

-  \  J  [ReOT^lMu11)  +  Im(Tm)Iffl(un)]  ds  (A-5) 

s 

After  multiplying  through  by  two,  this  can  be  written  as 

Re  (Tnum)  ds  -  Re  (Tmun)  ds  (A-6) 


for  this  to  be  true,  then 


or 


Re(Tnum)  =  Re(unTm) 
-ReClV1)  -  -Re  (u11!™) 


(A-7; 

(A-8) 


If  the  matrix  element  is  defined  as  um  then  from  Eq  (A-8) 


^nm  * 


(A-9) 
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Appendix  C 

Description  of  Computer  Program 


The  computer  program  in  this  appendix  is  for  the  analysis  of  the 
vibration  characteristics  of  an  open  cylindrical  cantilevered  shell. 

The  program  consists  of  the  executive  program.  Shell,  followed  by  two 
subroutines.  Fill  and  Poly.  The  purpose  of  the  program  is  to  calculate 
the  elements  of  the  system  matrix,  Q,  for  a  given  frequency,  u,  and  then 
calculate  the  matrix  determinant.  The  vanishing  of  the  determinant  es¬ 
tablishes  that  a  natural  frequency  of  the  system  has  been  found.  The 
program  has  not  been  written  to  Iterate  to  a  solution,  but  to  only  sweep 
across  a  range  of  frequencies  at  a  desired  step  size,  while  printing  out 
the  determinant  after  each  calculation.  Since  the  program  can  be  run 
interactively,  the  operator  can  choose  the  input  parameters  to  manually 
iterate  as  close  to  a  solution  as  desired.  A  description  and  listing  of 
the  program  and  subroutines  follow. 


Program  Shell 

Program  Shell  is  the  executive  routine  which  reads  in  the  required 
input  parameters,  and  then  after  calling  the  appropriate  subroutines, 
uses  a  system  subroutine  to  calculate  the  determinant  of  Q.  The  program 
does  this  a  specified  number  of  times  over  a  range  of  frequencies. 


Subroutine  Poly 

Subroutine  Poly  calculates  the  values  of  X^,  <(>^  and  and 

stores  them  in  array  C.  It  performs  these  calculations  for  n  ■  0  to 
?!v'  NTERM-1.  It  also  fills  the  array  CJ  with  the  complex  conjugate  of  the 


70 


values  in  array  C.  This  subroutine  utilizes  two  system  routines.  One 
is  needed  to  solve  an  eighth  order  equation  for  its  roots,  and  the 
other  to  solve  a  system  of  complex  equations  for  their  unknown  coeffi¬ 
cients  . 

Subroutine  Fill 

Subroutine  Fill  uses  the  C  and  CJ  arrays  to  calculate  the  elements 
of  the  Q  matrix.  The  coding  of  the  program  is  representative  of  Eqs 
(52)  and  (54)  in  the  text. 
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This  report  develops  a  variational  technique  for  the  analysis 
of  the  vibration  characteristics  of  an  open  cylindrical  cantilevered 
shell.  The  technique  is  developed  by  modifying  Reissner's  principle, 
which  normally  applies  to  static  problems,  through  the  use  of  Hamil¬ 
ton's  principle  so  that  it  applies  to  dynamic  problems.  The  varia¬ 
tional  technique  is  first  derived  in  general  for  an  elastic  system, 
and  then  specifically  tailored  to  an  open  cylindrical  cantilevered 
shell.  The  technique  is  implemented  by  first  finding  a  general 
solution  which  satisfies  the  equations  of  motion  for  a  cylindrical 
shell.  A  method  is  then  formulated  to  use  this  general  solution  to 
construct  a  set  of  trial  solution  functions.  With  the  variational 
method,  the  coefficients  to  this  trial  solution  function  are  then 
calculated  so  that  the  function  not  only  satisfies  the  equations  of 
motion,  but  also  the  boundary  conditions  around  the  four  edges  of  the 
shell.  A  computer  method  was  developed  to  perform  the  necessary  cal¬ 
culations  to  implement  the  variational  procedure,  but  preliminary 
results  have  shown  that  numerical  problems  must  be  eliminated  before 
accurate  results  can  be  expected. 

Experimental  data  for  an  open  cylindrical  cantilevered  shell  was 
also  collected  on  a  modal  analyzer.  The  results  are  presented  and 
discussed. 
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